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Abstract 

For the hopping dynamics in a one-dimensional model, containing energy and barrier disorder, we 
determine the linear and nonlinear response to an external field for arbitrary external frequencies. 
The calculation is performed in analytical terms. We systematically analyze the parameter space 
and find three different regimes, corresponding to qualitatively different frequency dependencies 

i 

of the nonlinear response. Two regimes agree with the results of recent conductivity experiments 
on inorganic ion conductors and ionic liquids, respectively. The ratio of the nonlinear and linear 
conductivity in the dc-regime can be explicitly expressed in terms of the disorder parameters. As 
a generic feature the nonlinear conductivity displays a minimum as a function of frequency which 

00 

can be identified with forward-backward dynamics in a double-well potential. The magnitude and 
sign of the nonlinear conductivity around the minimum is a measure of the disorder, inherent in 
this model. Surprisingly, the frequency of the minimum is hardly influenced by the disorder. 
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I. INTRODUCTION 



Nonlinear ion transport in disordered systems, like glasses and liquids, is of interest for 
mainly two reasons: (i) In electrochemical devices, the integration of thin-film electrolytes 
reduces the overall electrical resistance. In thin films, even small voltage drops may lead 
to high electric fields and to a field-dependent ionic conductivity, (ii) From a basic science 
point of view, the study of nonlinear effects yields additional information about ion transport 
mechanisms. 

When high dc electric fields E^c are applied to isotropic ionic conductors, the field de- 
pendence of the nonlinear current density jdc(E dc ) can be described in a first approximation 
by m\ : 

Here, q denotes the ionic charge and a app the 'apparent jump distance'. 

A experimental method for measuring nonlinear conductivities is the application of large 
ac electric fields E(t) = E cos(ut). In the dc-limit E can be identified with E dc - Non- 
linear conductivity then leads to the presence of higher order harmonic terms for the re- 
sulting current. For isotropic systems nonlinear effects are reflected in terms proportional 
to cos(3o;t), sin(3o;t) and higher harmonics. Taking into account only those terms which 
represent conductivity in phase with the electric field one can generally write [7]: 

j(t) = (a[(u) + 0{El)) E cos(ut) + \ (a' 3 (u) + O(E 2 )) E 3 cos(3ut). (2) 

In general cr^(cu) is frequency dependent for disordered systems. For the dc-limit the 
validity of EqjT] translates into the following expressions: cr[(uj — > 0) = joq a a,pp/ {2k B T) and 
a' 3 (co -> 0) = Jo(qa app /(2k B T)) 3 /6. 

Via the relation 

a[(uj ->■ 0) 6 \2k B TJ U 
it is possible to determine the apparent hopping distance a app from knowledge of cr' 3 (u — > 0) 
and <y[(uj —> 0). 

Typically, measured values for a app range between 1.5 nm and 3 nm p~H6]. These 
values are much larger than the typical hopping distance a hopp of ion transport which, e.g., 
for alkali silicate systems is around to 0.25 nm [8]. It has been shown theoretically that, in a 
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disordered potential landscape, the value of a app cannot be related to aho PP in a simple way 
[7J . Only in a regular potential, a app is identical to the hopping distance ah opp of the ions. 

Experimentally, nonlinear ac measurements have been carried out on alkali-ion conduct- 
ing glasses as well as on the ionic liquid HMIM-TFSI [71 19HTB"]. In both cases, the nonlinear 
conductivity a' 3 was positive in the dc regime and gradually decreased with increasing fre- 
quency. In the case of the alkali ion conducting glasses, a 3 became negative in the dispersive 
regime. Furthermore, the onset frequency for the linear and nonlinear conductivity, from 
which on deviations from the dc-regime become relevant, are nearly identical. In contrast, 
a' 3 remained positive for the ionic liquid over the entire frequency range. Furthermore it 
turned out that the frequency-dependence already starts at lower frequencies for the nonlin- 
ear conductivity [13]. Finally we mention that from general arguments it follows (see, e.g. 
[7J) that one expects a positive high-frequency regime. This has been explicitly shown, e.g., 
for simulations of a one-dimensional hopping model [TJ. 

The theoretical understanding of the ion dynamics is complex due to the strong interac- 
tion among the mobile entities [T1HT7] . Several groups have studied single-particle hopping 
motion in a discrete disordered energy landscape to learn about the linear conductivity; 
see, e.g., [T8W2T] . Recently, this approach has found a numerical justification since to a 
good approximation the ion dynamics can be mapped on a single-particle vacancy dynamics 
between distinct sites [22J. 

For ID models it is possible to analytically calculate the linear and nonlinear dc-current 
PHI 1231 [21] • Interestingly, for periodic boundary conditions the conductivity displays non- 
analytic behavior at zero field in the thermodynamic limit [7] El]. From studying a disordered 
discrete 3D-energy landscape via computer simulations it turned out that for a Gaussian 
distribution of site energies the experimental relation a app ^> ah opp can be recovered [25] . 
Surprisingly, for box- type distribution <j' 3 (uj — > 0) is even negative, so that according to Eqj3] 
it is not possible to get a real value for ah opp - Thus, a priori the observation a app ^> aho PP is 
non-trivial. 

To understand the frequency-dependence of cr' 3 (u), the dynamics within a double- well 
potential were studied in [26J. Except for the positive high-frequency plateau, c 3 (uj) is 
negative throughout the remaining frequency range. In general it is not evident whether 
the restriction to a double-well potential is sufficient to explain the frequency-dependence 
of a' 3 (u). 
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Here we present a periodic ID model which, on the one hand, displays typical behavior 
of disordered systems but, on the other hand, is simple enough so that the frequency depen- 
dence of the linear and the nonlinear conductivity can be calculated analytically. Discussion 
of this model will broaden the understanding of the experimental observations described 
above. This model contains three parameters, reflecting the typical barrier heights, the site 
disorder and the barrier disorder. Two questions are of particular relevance: (1) How do 
these parameters influence the qualitative frequency-dependence of the nonlinear conduc- 
tivity? (2) To which degree can the theoretical description of the minimum in the nonlinear 
conductivity be reduced to the dynamics in local double-well potentials? 

The paper is organized as follows. In Sect. II we introduce the model. The analytical 
results can be found in Sect. III. They will be discussed in Sect. IV. We conclude in Sect. V. 

II. MODEL 




FIG. 1: Sketch of the hopping model, used in this analysis. 

As a model we consider a periodic potential with two different wells. It is characterized 
by the three energy values Vi, V 2 and V3; see Figjl} In what follows we always consider 
V\ < V 2 . A particle performs activated hopping between adjacent sites. We introduce the 
probabilities 2/1,2 (£) that the particle is in the left or the right well, respectively. In particular 
one has yi(t) + 2/2 (i) = 1. The activation barriers are modulated by the external potential 
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Vextif) = ecosut. 

In case of a regular potential, i.e. V\ = V 2 and V 3 = this simple model yields 

j dc oc sinh (^j . (4) 
Via comparison with Eqjl] we can relate e to the experimental parameters, i.e. 

£ = qa^ppEo/2. (5) 

As mentioned above, for this regular potential one has a app = a hopp . For reasons of 
simplicity we will express all our results in terms of e. Furthermore, we choose fcgT= 1, i.e. 
we express all energies in terms of k^T . 

There are two possible transitions from well 1 to 2 and vice versa. This yields for the 
transition rate from well 1 to well 2 

= ri_> 2 + r 2 ^i (6) 

with 

IW = r exp(-(V r 1 + V 3 ) - e cos(wt)) (7) 

and 

r^ 2 = T exp(- V 2 + e cos(w*)) (8) 

For simplification we choose the kinetic prefactor r to be unity. Introducing the notation 
Ti := e * one can therefore write 

1 12 — i ie + 1 ii 3 e . {t)) 

In analogy one obtains 

r 21 = r 2 e- £COSW< + r 2 r 3 e ecoswt . (10) 

Now one can write the rate equation for yi(t) as 



vx{t) = -r 12 {t) yi {t) + r 21 (t)(i - yi {t)) 

= - yi (t) ((Ti + r 2 r 3 )e £COS ^ + (r 2 + r 1 r 3 )e" £COS ^) 
+ r 2e - £COSWt + r 2 r 3 e ecos ^ 



(ii) 



(12) 
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To simplify this rate equation we introduce the equilibrium population of yi, denoted as 
Pi, given by the Boltzmann factor (AV = V<i — V\) 

e -AV/2 
11 '~ gAV/2 _|_ e -AV/2 

(13) 

and define u\{t) as the difference to the equilibrium population, i.e. 

ui(t) := Vl (t) - Pl . (14) 
Then a straightforward calculation [27] yields from Eqjl~2| 

u^t) = - Ul {t)(r a e £COSUJt + r b e- £COSa ") + r^- e £COS ^ + e- £COSW '). (15) 
Here we used the abbreviations T a := Fx + r 2 r 3 and Tj, := T 2 + TiT 3 as well as 

_ e -(l/2)(Vi+V2)-V 3 _|_ e -(l/2)(Vi+V 2 ) 
e AV/2 _|_ e -AV/2 

(16) 

Note that after a transformation V\ — > V\ — c, V% — > Vi — c and V3 — > V3 with c < V\ and 
c < V2 the dynamics is the same except for a trivial scaling of the hopping rate by exp(c). 
Thus, in the subsequent analysis we can set V\ = 0, i.e. F\ — 1. Then, V2 can be interpreted 
as the asymmetry, V3 as the barrier disorder, Eqnjl6] reduces to 

r r 2 (i-r 3 ) 



1 • v, ' 

and EqnfT3~1 to 



(17) 



» ^ rk- (18) 



III. ANALYTICAL SOLUTION 



A. Non-equilibrium population 



The general solution of can be written as a sum over the different harmonics cos(no;t) 
and sm(ncut), i.e. 

00 00 
Ui(t) = a n (e) cos(nut) + '^^b n (e) sm(nut). (19) 

n=0 n=l 

Formally, the individual a n (e) and b n (e) can be written as a Taylor series in e. One can 
easily show that the lowest-order terms of a n (e) and b n (e), respectively, are proportional to 
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To calculate the current which is in phase with the external field only the cosine-terms, 
i.e. the a n (e), are of relevance. In what follows we define 

a n := lim{a n {e)/e n ). (20) 

e— >0 

Since for small fields a n oc e n , the a n are independent of e and characterize the nonlinear 
dynamics in the experimentally relevant case of weak nonlinear effects. Similarly, one can 
define the (3 n as new coefficients for the sine terms. 

After inserting Eqjl9] into EqJT5| and expanding the exponential terms in terms of the 
different harmonics (cos(no;t), sin(nu;t), obtained after application of appropriate addition 
theorems) one can set all terms, belonging to the same harmonics and the same exponent 
in e, to zero. This yields a system of linear equations in the a n and (3 n which can be solved 
iteratively, starting with n = 1. After a tedious but straightforward calculation one ends up 
with 



«1 = ^X~2' W 

to 2 + % 



and 



-T(2u 2 - 7 2 ) 

" 2 ~ 1 - (u; 2 + 7 2)(4u; 2 + 72) (22) 

7 + r( 7 !-5a; 2 ) ^{iW-^) 

Qs 6(9w 2 + 7 2)(a; 2 + 72) +7 -2(w 2 + 7 2)(4w 2 + 72)(9c l ; 2 + 72)' 1 ' 



Here we use the abbreviations 



7+: =r a + r 6 = (i + r 2 )(i + r 3 ) (24) 
7 _:=r -r 6 = (i-r 2 )(i-r 3 ). (25) 



B. Calculation of the conductivity 

The dynamics of a charged particle in the periodic potential can be characterized by a 
current. Based on the time-dependent populations as determined in the previous section 
one can generally write for the current 

j(t) = rw^VW - r 2 e- £COS ^ 2 (t) + r 2 r 3 e ecos ^ 2 (t) - r 1 r 3 e- £COS ^i (t) (26) 

which can be rewritten with ?/ 2 (t) = 1 — yi(t) as 

j(t) = -r 2 e- £COSWt + r 2 r 3 e ecoswi + yi (t)(r c e £COSU}t + r>- £COS ^). (27) 
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Here we have used the abbreviations T c := I\ — and Td := T2 — T1T3. From now on 
we will choose again r\ = 1. 

For the subsequent calculations the exponential-terms have to be expanded with respect 
to e. Furthermore, for the population yi(t) the expression from Eq, 19 combined with Eqjl4| 



can be inserted. Then one needs to combine again all terms which scale with the same 
harmonics and the same exponent in e in analogy to the calculation of the a n . This results 
in an expression for the current which contains terms proportional to cos(no;t) or sin(no;t) 
for all integer n > 1, as well as the coefficients a n , f3 n . 

Macroscopic disordered systems behave isotropically, i.e. the reversal of an electric field 
leads to a reversal of the current. In contrast, the present model is anisotropic for V2 7^ 
0. To enable a direct comparison with the experimental situation, we therefore average 
over two opposite directions of the field. As a consequence only the odd harmonics (i.e. 
cos(ut), cos(3o;t), ...) and only terms which are uneven with respect to e remain. All other 
terms cancel. Thus, one can finally write 

j{t) = Si(e, u) cos(cot) + s 3 (e, uj) cos(3wt) + sine terms, (28) 

where the lowest order term of s n (e,u) is proportional to e n . In analogy to the discussion 
of the populations a n we define the conductivities via 

o- n (w) :=Um(s n (e,w)/e n ), (29) 

thereby capturing the effect of the external potential in lowest order. 

Via comparison with Eqj2]one can identify o\ = a[ and 03 = CX3/4. The experimentally 
relevant quantity a a pp/flhopp can thus be expressed as (also using Eqg 

a app _ 24a 3 (u = 0) 



ahopp (Xi(u; = 0) 

More generally we define 



(30) 



A ^ - (31) 
In what follows we give the results of the lengthy but straightforward calculations de- 
scribed above [27]. The in-phase terms read 

= ] { M r 3) + (1 + T 2 )(l - T 3 ) ai (u) (32) 

1 T J- 2 
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and 



= I |^-p r g + |ai(o;)(l + r 2 )(l-r3) + ia 2 (a;)(l-r 2 )(l + r3) 

+ a 3 M(i + r 2 )(i-r 3 ). (33) 



These equations are the key result of this work because all conclusions, discussed below, 
follow from them. 

Of particular interest are the limits at vanishing or infinite frequency. After some algebra 
one obtains 

*<" = > = ( i + £aW (34) 



cr 3 (u; = 0) 



^ = 0) A , r (i-r 2 ) 2 (i-r 3 ) 2 \ ,,,, 
-^^l 1 + 6 (i + r 2 ) 2 (i + r 3 )0' (35) 

a 1 (^oo) = 2 ^±^, (36) 

1 + 1 2 

and 

a 3 ( W -> oo) = -^-^ (37) 

Furthermore, it is also possible to express the whole frequency-dependence of the linear 
conductivity by a simple expression. It is given by 

a 1 {oj) = (7 1 {u = 0) + 2 , 4 2 F3 (38) 

7+ + w z 

Naturally, this expression is consistent with the two limiting cases, given above. 





FIG. 2: Comparison of the numerical solution of the rate equations (symbols) with the analytical 
expressions (lines) for F 2 = 2 and different values of V 3 . Left: cri(w), right: <t 3 (w). For labeling 
purposes in Fig(3]the three parameter specifications are denoted a, b, c from top to bottom. 
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Due to the complexity of the expressions we first check that straightforward numerical 
simulations of the corresponding rate equations and subsequent determination of the current 
via Fourier transformation agree with the analytical expressions. This is shown in Fig|2j 
Indeed, one can find an excellent agreement (also, for other values of V2 and V3 not shown 
here) which strongly supports the correctness of our algebraic calculations. 



IV. DISCUSSION 

First, we discuss A(u) as introduced in Eqj3l| For the special case of an ordered potential 
(i.e. V2 — V3 — O), one obtains A(u) = 1. Generally, in the high-frequency limit all values a n 
disappear, i.e. the populations of the two wells do not differ from the Boltzmann distribution. 
As a consequence, A(u — > 00) is unity. One can also see from Eqsj21~|-|2~3|that for the specific 
case T 3 = 1 all a n disappear for all finite frequencies. As a consequence, both the linear and 
the nonlinear conductivity are frequency-independent and thus trivially display the same 
ratio for all frequencies. This limit corresponds to the case of vanishing barrier disorder, i.e. 
V3 = 0. A prototype model for this scenario is the trap model [28] . 

Of particular interest is the zero-frequency limit of A(oS). Based on Eqs,32 and 33 one 
obtains 

Thus, for V2 = or V3 = , equivalently T 2 = 1 or T 3 = 1, the proportion between 
the nonlinear and linear part of the conductivity becomes minimal. For large values of 
the asymmetry and the barrier disorder, A(u = 0) can approach values as large as 7. 
Thus, for this model any disorder generally leads to A(u = 0) ^> 1 in agreement with the 
experimental observation. Interestingly, for the dc-behavior the impact of the asymmetry is 
the same as the impact of the barrier disorder. Naturally, if expressing these results in terms 
of temperature one would expect an increase of A(u = 0) with decreasing temperature. 

It turns out that in the whole parameter space 173(00) displays one minimum as a function 
of frequency. However, as shown in Fig|2] , two other properties , i.e. existence of negative 
values and sign of initial slope of u 3 (u) strongly depend on the specific values of V 2 and 
V3 as shown by variation of V3 for fixed V2 (here: V2—2) . First, for V3 = 1.2 the nonlinear 
conductivity does not acquire negative values whereas for V 3 > 1.8 negative values are 
observed. Second, the low-frequency slope of 03(0;) is negative for V3 < 1.8 and positive else. 
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FIG. 3: Representation of the different regimes I, II, and III (definition, see main text) concerning 
the qualitative behavior of the frequency dependence of a^icj). The letters a,b,c mark the values 
V2 and V3 that are used for the three curves in Figj2j The numbers 1-4 refer to the different limit 
cases mentioned in the text. 

This shows that, for certain V2, all possibilities can be realized through variation of V3. Via 
a numerical analysis we have analyzed each parameter pair (V2, V3) with respect to these two 
properties. We have identified three different regimes: (I) Positive initial slope and presence 
of a negative frequency regime; (II) Negative initial slope and presence of negative frequency 
regime; (III) Negative initial slope and only positive values. The remaining possible option 
does not occur. The experimental data for ion conductors correspond to regime (II), those 
for ionic liquids to regime (I). A systematic identification of all three regimes can be found 
in Fig|3j 

Interestingly, the presence of a negative regime hardly depends on the specific value of 
the asymmetry V2 but rather on the difference of the barrier heights V3. In contrast, regime 
I is approximately symmetric with respect to the values of V2 and V3. 

Some of the limiting values of this phase diagram can be predicted even analytically. Due 
to the complexities of the resulting algebraic equations we have employed an algebra program 
(Mathematica). The following results were obtained: (1) The transition between regime II 
and III for V2 = occurs for V 3 = ln(|(29 + 4\/7)) « 1.48. (2) The transition between 
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FIG. 4: Left: 03(0;) for fixed V2 = 1 and for different V3 G [2,6], increasing from top to bottom in 
steps of 0.5. Right: The same as before but with a shifted y-axis. 

regimes I and II for large V3, i.e. for vanishing slope of 03(0;), occurs for V2 = In 3 ~ 1.10. (3) 
The transition between both regimes for large V2 occurs for V3 = ln((73+10v^46)/27) ~ 1.65. 
(4) The transition between regimes II and III for large V2 occurs for V3 ~ 1.62. This last 
value emerges from solving a polynomial equation of 12th order and was determined only 
numerically. 

In the limit V3 — > 00, i.e. 1^3 = , the periodic potential transforms into a double-well 
potential DWP with vanishing dc-conductivity, so cri(u; = 0) = 0-3(0; = 0) = 0. In this 
case the system can either be described by regime I or by regime II. In both cases 03 (u) 
becomes negative for intermediate frequencies. Most importantly, upon variation of (mainly) 
V3 the nonlinearity does acquire negative values (large V3 , regime I or II) or only positive 
values (small V3 , regime III). Thus, variation of the sign of <Js(u) at its minimum does not 
necessarily reflect different physical mechanisms. In this sense the different experimental 
results for inorganic ion conductors and ionic liquids mentioned in the Introduction may 
just reflect a different degree of disorder. 

One may wonder whether the occurrence of the minimum in 03(0;) can be related to 
forward-backward motion in a local double-well potential , which would not contribute to 
an overall conductivity. If this were the case, the limit 1/3—7-00 would already contain the 
relevant information about this frequency regime. In order to elucidate this aspect , we have 
calculated 03(0;) for a fixed value of V% and various values of V3. In order to stay in the 
experimentally relevant regime II we have chosen V2 = 1. The resulting graphs are shown 
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FIG. 5: The characteristic frequencies uj\ and ui^ as a function of V3 for fixed V2 = 1. The thin 
solid line corresponds to the approximation Eq|41| 

in Figji] (left). Interestingly, the high-frequency plateau is nearly unchanged by V3 and thus 
indeed just reflects the properties of the local DWP, governed by the asymmetry V2. This 
can be directly read off from the analytical expression EqfSTj However, also the frequency 
regime left of the minimum w«l displays some specific properties. If one subtracts cr 3 (0), 
the dc-value of the nonlinear conductivity, one finds a very similar behavior for different 
values of V3 in this low- frequency regime to < 1; see Figji] (right). Thus, in this frequency 
regime the nonlinear conductivity is basically a sum of the dc-value and the contribution 
of the double-well potential, described by a 3 (u,V 3 = 00). It turns out , though, that this 
simple superposition principle becomes worse in regime I or for V% ~ 0. 

For a closer understanding of the nonlinear conductivity it may be of interest to compare 
its characteristic frequencies with those of the linear conductivity. We start by compar- 
ing the two onset-frequencies uj* and u> 3 which characterize the frequencies from which on 
deviations from the dc-regime are relevant, i.e. where the conductivity increases or de- 
creases, respectively. Here we consider an increase and decrease of 10%, respectively, i.e. 
&x(u%)/<Ti{u = 0) = 1.1 and <j$(u%) / <j$(u = 0) = 0.9. The qualitative behavior does not 



depends on this specific choice of 10%. w* can be calculated analytically from Eq.38 One 
obtains 

= 7+ 1 _ p (40) 
which for large V3 (i.e. T3 <C 1) can be approximated as 

ul « (1 + r 2 ) y/0W 3 oc exp(-l/ 3 /2) (41) 
13 
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FIG. 6: The characteristic frequency oo™ m as function of Vi and V3. 

For the numerical analysis we restrict ourselves again to V% — 1 in order to avoid the 
regime III for better comparison with the experimental situation. We have varied V3 over a 
broad range, encompassing the regimes I and II. As shown in Fig|5] there exists a significant 
dependence of u\ on V3. The approximation Eq|4~Tj i.e. the scaling with exp(— V3/2), roughly 
works for V3 > 3. 

Due to the complexity of 03 (cu) one cannot write down a simple analytical expression for 
Therefore we restrict ourselves to the numerical determination of bj\] see Fig|5j Two 
important observations can be made: (1) The V^-dependences of u\ and u\ are identical for 
large V3 (regime II). Both frequencies only differ by a constant factor close to unity. (2) For 
small V3, i.e. in regime I, one observes Cw*. We would like to stress that this is fully 
compatible with the experimental observations as reviewed in the Introduction (inorganic 
ion conductors: regime II; ionic liquids: regime I). 

Another characteristic feature of cr^(u) is its minimum, occurring at a frequency oj™ m . 
The data in Fig|2] already suggest that u™ m has no strong dependence on V3. We have 
analyzed its dependence on V2 and V3. In Fig|6] we show more systematically, how uj™ %n 
depends on V2 and V3. Obviously, the dependence is very minor; less than a factor of two if 
taking into account all possible parameter pairs (V2 > 1, V3). 

This has several interesting implications: (1) Whereas variation of V2 and V3 does not 
modify oj™ in si gnificantly, this is not true for variation of V\. As discussed above consid- 
eration of the Vj-dependence gives rise to a trivial factor exp(— Vi) for all rates and all 
characteristic frequencies such as o;™ m . This exponential dependence on V\ thus has a sig- 



14 



nificant impact on the frequency scales of the conductivity. Thus, the minimum frequency 
of the nonlinear conductivity is mainly sensitive to the typical barrier height between ad- 
jacent minima rather than the asymmetry or the barrier disorder. (2) Comparison of the 
linear and the nonlinear conductivity in FigfJ] leads to the assumption that the appropriate 
counterpart to u™ m is the frequency where the slope of <J\(u) is maximum in the double- 



logarithmic-representation. We denote this frequency by co™ ax . Starting from Eq.38 it can 
be calculated analytically. One obtains 

«!- - y+ (^) V4 - a + r.)d + r.) (^) " 4 ■ (42) 

Since the dependence on T 2 is only via the factor (1 + r 2 ), the dependence on T 2 vanishes 
for large V 2 in agreement with u™ m . In contrast to w™ m the frequency uj™ ax exponentially 
depends on V3 for large values of V3 via exp(— V3/4). However, due to the factor (1/4) in 
the exponent, the V 3 dependence is significantly weaker than for us\ oc exp(— V3/2). Thus, if 
at all, the minimum of the nonlinear conductivity is related to the region of maximum slope 
of the linear conductivity as evaluated in a double-logarithmic representation. 



V. CONCLUSION 



The analyzed ID hopping model may be considered as a minimum model, which captures 
the non-trivial frequency dependence of the nonlinear conductivity. Despite its simplicity, 
the algebraic calculations are quite lengthy and some results required the help of an algebra 
software. A key result was the presence of a minimum of the nonlinear conductivity for 
all parameters except for the trivial ordered case. The scaling properties of the different 
characteristic frequencies of the linear and nonlinear response allow one to see how the 
disorder influences the shape of the frequency-dependent conductivities. In comparison 
with the experimental results on inorganic ion conductors and ionic liquids it seems that the 
disorder effects are more pronounced in the first case. 

Of course, this model mimics the true experimental system only in a very simple way. 
First, one has to consider that in reality one has a 3D rather than a ID energy landscape. 
For example it turns out that for a random barrier landscape, here corresponding to V2 = 
one has A(u = 0) = 1 in the ID system as discussed above but A{u = 0) < for the 3D 
system [26] . Second, in reality a disordered system is likely to be described by a distribution 
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of energy parameters rather than by well-defined values. Of course, a natural next step 
would be to average the present results over appropriately chosen parameter distributions 
as performed, e.g., in [2o] . 

Despite these limitations one may expect that important general features of the model 
may also hold for more complex systems. In particular the observation that the frequency 
dependence close to the minimum of the nonlinear conductivity in the experimentally rele- 
vant regime is to a large extent related to the dynamics in a local DWP may justify the use 
of the present model for a qualitative understanding of experimental systems. 

Acknowledgements 

We gratefully acknowledge the support by the DFG (FOR 1394) as well as very helpful 
discussions with M. Lowe. 



[1] J. Vermeer, Physica 22, 1257 (1956). 

[2] L. Zagar, E. Papanilolau, Glastechn. Ber. 42, 37 (1969). 

[3] J. P. Lacharme, J. O. Isard, J. Non-Cryst. Solids 27, 381 (1978). 

[4] J. M. Hyde, M. Tomozawa, Phys. Chem. Glasses 27, 147 (1986). 

[5] J. L. Barton, J. Non-Cryst. Solids 203, 280 (1996). 

[6] J. O. Isard, J. Non-Cryst. Solids 202, 137 (1996). 

[7] B. Roling, S. Murugavel, A. Heuer, L. Liihning, R. Friedrich, S. Rothel, Phys. Chem. Chem. 
Phys. 10, 4211 (2008). 

[8] R.D. Banhatti, A. Heuer, Phys. Chem. Chem. Phys. 3, 5104 (2001). 

[9] S. Murugavel, B. Roling, J. Non-Cryst. Solids 351, 2819 (2005). 

[10] A. Heuer, S. Murugavel, B. Roling, Phys. Rev. B. 72, 174304 (2005). 

[11] H. Staesche, B. Roling, Phys. Rev. B 82, 134202 (2010). 

[12] H. Staesche, B. Roling, Z. Phys. Chem. 224, 1655 (2010). 

[13] F. Feiten, B. Roling, Solid State Ionics 226, 59 (2012). 

[14] M.D. Ingram, Phys. Chem. Glasses 28, 215 (1987). 

[15] J.C: Dyre, P. Maass, B. Roling, D.L. Sidebottom, Rep. Prog. Phys., 72, 046501 (2009). 



16 



[16] P. Maass, J. Non-Cryst. Solids 255, 35 (1999). 

[17] R.D. Banhatti, K. Funke, Solid State Ionics 177, 1551 (2006). 

[18] A.G. Hunt, J. Phys. Cond. Mat. 3, 7831 (1991). 

[19] S.D: Baranovskii, H. Cordes, J. Chem. Phys. Ill, 7546 (1999). 

[20] J.C. Dyre, T.B. Schr0der, Rev. Mod. Phys. 72, 873 (2000). 

[21] T.B. Schr0der, J.C. Byie,Phys. Rev. Lett. 101, 025901 (2008). 

[22] H. Lammert, A. Heuer, Phys. Rev. Lett. 104, 125901 (2010). 

[23] K.W. Kehr, K. Mussawisade, T. Wichmann, W. Dieterich, Phys. Rev. E 56,R2351 (1997). 

[24] M. Eimax, M. Koerner, P. Maass, A. Nitzan, Phys. Chem. Chem. Phys. 12, 645 (2008). 

[25] S. Rdthel, R. Friedrich, L. Luhning, and A. Heuer, Z. Phys. Chemie 224, 1855 (2010). 

[26] B. Roling, J. Chem. Phys. 117, 1320 (2002). 

[27] C. Mattner, BSc Thesis, University of Minister (2012). 

[28] J. P. Bouchaud, A. Comtet, C. Monthus, J. Phys. I France, 5, 1521 (1995). 



17 



